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Stability of free-floating ship 


Part II 


Maciej Pawtowski 
Gdansk University of Technology 


ABSTRACT 


This is the second part of the paper published in Polish Maritime Research no. 2/2005, 
dealing with the calculation problem of righting arms of the free-floating ship, i.e. longitu- 
dinally balanced at any heel angle. In such case the righting arms are ambiguous as they 
depend on a way the heeling moment acts. Two cases were considered : when the heeling 
moment is parallel to the ship plane of symmetry, and the case when it performs the least 
work, i.e. when the moment is parallel to the main axis of ship waterline. It was demon- 
strated that angular translations (heel and trim) are then the Euler angles associated with 


a relevant reference axis. Some cases of the incorrect defining and using of those angles in today design 
practice were indicated. The most important features of the curve of righting arms of free-floating ship 


were demonstrated. 


Keywords : stability, free floating ship 


FEATURES OF RIGHTING ARM CURVE 


Knowing metacentric radii of a free-floating ship one can 
easily find the remaining characteristics of its curve of righting 
arms. They are analogous to those known from the classical 


theory of ships. 
Hence, the metacentric height A is equal to : 
h = Yap GZ = BM - BZ (25) 
where : 
BM - metacentric height according to Eq. (12) 
BZ=-rn - height of ship centre of gravity over its centre of 


buoyancy (Fig. 2) 

radius-vector of ship centre of buoyancy relative 
to its centre of gravity, and 

versor normal to ship waterline according to 
Eq.(6) or (20). 


Eq. (25) can be obtained immediately by considering the 
buoyancy action line in the rotation plane (Fig.6) for the heel 
angle increased by dọ. The metacentric height can be also ob- 
tained by differentiating the righting arm GZ, acc. to Eq.(10), 
respective to the heel angle in the ship-fixed reference system. 

The derivative is given by the formula : 


GZ' = e'-(rxn) + e-(r'xn) + e-(rxn') = BM + r-n 


- 
II 
Q 
w 


5 
1 


where the sign [' ] stands for differentiating respective to 0. 
It can be demonstrated that the first term equals zero, the 
second gives the metacentric height BM, and the third is equal 
to r:n. Hence the formula is identical with Eq. (25). 


Fig. 6. Ship rotation plane 


Work done by the heeling moment M 
is expressed by the formula : 


o ty) 
L=|Mdo=A[GZdo=Al, (26) 
0 0 
where : 


A - ship buoyancy 
lą - dynamic arm. 
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It can be observed that the area under the curve of righting 
arms, called the dynamic arm, is proportional to the work of 
heeling moment. 

Considering rotation of the rotation plane (Fig.6) by the an- 
gle do, one can easily demonstrate that the differential GZ do = 
= d(BZ) is an increment of the segment BZ due to vertical 
translation of the point Z, as the buoyancy centre B translates 
horizontally relative to ship's waterline. Hence the known 
formula for the dynamic arm lq = BZ — a, is obtained, where 
a = BoG - the height of the gravity centre G over the buoyancy 
centre of ship in the upright position (i.e. for @ = 0). The equa- 
tion has simple physical interpretation — the dynamic arm is 
equal to the vertical increment of the distance between the cen- 
tre of gravity and centre of buoyancy. It can be useful in chec- 
king the calculation accuracy of the righting arm curve. 

Worth mentioning that the curve of righting arms of free- 
-floating ship complies with the theorem of minimum potential 
energy, i.e. ship's heel (understood as a rotation of the plane of 
rotation) by a given angle demands work done to be minimum. 
This is an important feature of the curve. To prove it let's obser- 
ve that to incline the ship from its position of longitudinal 
equilibrium is not possible without applying a trimming mo- 
ment and doing an additional work that could increase its po- 
tential energy. Hence it results that righting arms of a free-flo- 
ating ship are at the most equal to or smaller than those of 
a fixed-trim ship, that is clearly illustrated in Fig.1. If it is not 
the case, it means that some errors appear in the calculation 
algorithm. 


Cross-curves of stability 


The hull form arm, i.e. the arm of buoyancy force relative 
to the initial location of the centre of buoyancy, shown in Fig.6, 
is given as follows : 


lp = GZ + a:sinod (27) 


For a free-floating ship which changes its trim during he- 
eling, the hull form arms depend on the height of centre of 
gravity over BP. Because of trim changes the idea of cross- 
-curves of stability does not apply, strictly speeking, to a free- 
-floating ship. Let's observe that z - axis fixed to the ship does 
not lie in the rotation plane which is inclined to it by the angle 
O. Projection of the axis z to the rotation plane agrees with the 
axis BoZ (Fig.6)" . Therefore, by shifting the centre of gravity 
along Oz-axis by the quantity Azg it moves before the rotation 
plane by the distance Ax? = Azgsin@. As a result the ship beco- 
mes unbalanced and it must trim in the vertical plane by the 
angle doz = Ax2/BM_,, where BM; is the longitudinal meta- 
centric radius at a given heel angle p. The ship can be balanced 
only by changing its trim, without changing its heel angle. In 
this case the relationship between both trim angles is the same 
as in Eq. (23) but without minus sign , i.e. doz = dOcos. Hence 
the trim correction is d@ = da2/cosd. From Fig.6 it can be sta- 
ted that the new righting arm is : 


GZ, = GZ — Azg-cos8-sind + (D/V) daz (28) 


where : D - moment of deviation 
of ship waterline in En - coordinate system (Fig.4). 


For normally occurring trim values the function cos can be 
omitted as practically being equal to 1. Two first terms of Eq. (28) 
are the same as for fixed-trim ship. The last term, (D/V) doz , 
denoted now by Al, accounts for the effect of trim on the righ- 
ting arm curve of a free-floating ship. Taking into account that 
daz = Ax7/BM_, one obtains : 


Al = Azg:sin0-tgx (29) 
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Hence it can be observed that it is possible to find the righ- 
ting arm curve for a new location of the centre of gravity by 
means of cross - curves of stability for free-floating ship provi- 
ded additional information on the run of the angles O and % in 
function of the heel angle ọ is available , which makes it possi- 
ble to calculate the correction Al. Let's observe that the correc- 
tion Al disappears not only when 6 = 0, as it could be expected, 
but also when % = 0, i.e. when the ship's floatation axis is paral- 
lel to the rotation axis. Moreover from Eq. (29) it results that 
the quotient Al/Azg is not dependent on a value of changes of 
the height of centre of gravity over BP; therefore in order to 
calculate the correction Al it is enough to know a course of the 
quotient Al/Azg = sinðtgy in function of the heel angle . With 
the use of the quotient, Eq. (28) for the new curve of righting 
arms obtains the following form : 


GZ) = GZ - Azg (sinb — AVAzg) 830) 


To improve calculation accuracy it is advisable that cross- 
-curves of stability are calculated in the form of the righting 
arm curves GZ for a selected, typical location of ship's centre 
of gravity. In such case the correction Al is small, hence it can 
be often neglected. Cross-curves of stability are usually pre- 
sented in the form of the diagram : GZ = GZ (V, = const). 

The diagram showing the differential quotient Al/Azg = 
= sin6tgx should be presented in the similar way, in function 
of the buoyancy V, for a fixed 9. 


OTHER DEGREES OF FREEDOM 


In ship hydrostatics in order to describe ship inclination 
relative to water, are used two mutually indpendent angles 
which can be interpreted as degrees of freedom, e.g. the analy- 
tical angles @ and O which describe orientation of waterline 
relative to ship. By using them all other angles between vario- 
us planes and axes can be expressed. 

However the analytical angles are not convenient to descri- 
be angular translations of ship's hull. To this end three natural 
angles (in dynamics called the Euler angles) are used. In static 
problems two angles are sufficient to describe angular transla- 
tions as any third angle is not necessary due to lack of yaw 
which means that nodal line orientation is constant in the sea 
plane. The above discussed angles and 0 associated with Oy- 
-axis and PS (playing the role of plane of nodes) exemplify the 
Euler angles. As there are two other axes, one can choose de- 
grees of freedom associated with the remaining axes. 

For instance, in NAPA software as well as many other com- 
puter programs’ the heel angle is defined as the slope angle of 
trace of water in stations, relative to Oy - axis, hence it is the 
angle @, whereas the trim angle is rotation of the ship around 
this trace, denoted by Ox . Hence they are Euler angles associa- 
ted with x-axis, commonly used in ship dynamics. The nodal 
line is the trace of water on midship section, fixed in space. 
The angle Ox is equal to that of slope of Ox-axis relative to 
water- level. As cos (90° + Ox) =i-n, then : sinx = cosatg6, or 
simpler : 

(31) 


The versor n can be directly expressed by means of the Eu- 
ler angles. Dividing sinx by tg@, one obtains : cosOx = cos&/ 
/cos@, hence : nz = cosa = cos0xcos@. By substituting them in 
Eq.(6), the components of the vector n take the following form : 


(32) 
As ny = - sing, hence : sino = cos@xsing, from which it 


results that : ¢<@, which has been already obtained from Eq.(7); 
both angles are equal to each other in case of no trim. From the 


tg0, = cos@-tgO 


n = (-sin®, , -cosOx-sin@ , cosOx:cos@) 


formula for cos it results that : @ < a, hence : d< a. The angle 
@ =q, in case of no trim or when the heel angle equals 90°. 


In the case of the reference axis Oz, 
the Euler angles are : 


> the inclination angle by which the axis z is deflected from 
the vertical, given by Eq.(5); it is identical with the slope 
angle of BP relative to waterline 
> as well as the trim angle y, i.e. the ship rotation angle aro- 
und Oz - axis. 
The edge of intersection of BP and water level forms 
a nodal line as its orientation is fixed is space. Both the angles 
and the nodal line is shown in Fig.3. The angle w is equal to the 
slope angle of trace of water in the plane z = const (1.e. nodal 
line) relative to Ox-axis . 


This way one obtains as follows : 
tgy = - tg0/tg@ 
From Eqs.(5) and (33) it results that tga can be taken as a vec- 
tor composed of tg@ and - tg0. Hence : tg@ = tgacosy, and : 


tg = - tgasiny. By substituting them in Eq.(6) the compo- 
nents of the vector n take the following form : 


(33) 


n = (siny-sina , - cosy-sina , cosa) 64 


As : ny = - sing, hence : sin = cosysina., from which it 
results that: ¢ < a; both the angles are equal to each other only 
in case of no trim. 

It is not indifferent which reference axis (and the Euler an- 
gles associated with it) is chosen to calculate righting arm cur- 
ve because location of the rotation plane in which the ship has 
to be balanced depends on the axis chosen. During balancing 
the ship, a change of trim angle does not influence the rotation 
angle. 

In the case of Ox reference axis the rotation plane is a ver- 
tical section perpendicular to water-level and parallel to water- 
-trace on frame sections (Fig.7); in the case of Oy-axis — it is 
perpendicular to water-trace on PS (Fig.2), and in the case of 
Oz-axis — perpendicular to water-trace on BP (Fig.8). In other 
words, in the case of the reference axis x the rotation plane is 
parallel to nodal line, and in the two latter cases — perpendicu- 
lar to nodal line. The three planes of rotation are not identical 
which makes that for a free-floating ship righting arms depend 
on the choice of reference axis. In all the cases the rotation 
plane is fixed in space, perpendicular to heeling moment vec- 
tor, vertical, and crossing ship's centre of gravity. The perma- 
nent orientation in space is associated with the permanent di- 


Fig. 7. Rotation plane of free-floating ship which trims 
around trace of water on midship section 


rection of nodal line in sea-level, due to the lack of yaw. Du- 
ring heeling the buoyancy centre moves along the curve of 
buoyancy centres located on the rotation plane whose features 
have been already discussed. For this reason the righting arm 
curve is a function of the rotation angle of rotation plane, fur- 
ther marked 1. In other words, the heel angle of the ship is 
identical with the rotation angle of the rotation plane. 


Fig. 8. Rotation plane of free-floating ship which trims around Oz- axis 


In the case of Ox-axis the rotation angle n cannot be iter- 
preted simply in geometrical terms. It is also difficult to realize 
such heeling model physically, that can be observed in Fig.7. 

As a result the reference axis is not located in the rotation 
plane. The rotation angle and heel angle are mutually related 
by, the formula : dn = d@cos6x, hence the rotation angle n = 
= Jcos0,d@, from which it results that : n < @, and that for ọ = 
= 90° the rotation angle of rotation plane is a little smaller than 
90° (n < 90°) . 

Generally, differences between the angles n and @ are mar- 
ginal as cos@, is practically equal to 1, nevertheless both the 
angles are not identical. For Oy-axis the rotation angle n is 
identical with the angle , and for Oz-axis — with the angle a. 
In two latter cases shown in Fig.2 and 8, the reference axis is 
parallel to the rotation plane and it rotates together with the 
plane. 

It can be mathematically demonstrated that the rotation pla- 
nes associated with particular reference axes are differently si- 
tuated relative to the ship. To this end it is sufficient to demon- 
strate that versors normal to the planes, called rotation axes, 
are different. 

When x-axis is a reference axis the vector normal to rota- 
tion plane is given by the formula : e = e2xn, where ez — versor 
of trace of water on frame sections. The versor e differs from 
the versor e; = (cos8, 0, sin) as for an arbitrarily inclined ship 
the water-traces on PS and midship section are not mutually 
normal. As e2 = (0, cos@, sin), hence the versor of rotation 
axis is as follows : 


(35) 


When Oz- axis is a reference one the versor e normal to 
rotation plane is parallel to the vector : kxn. Accounting for 
that the length of the vector |kxn| = sina, one obtains as fol- 
lows : 


e = (cosx , - sinOxsing , sin8,-cos@) 


(36) 


which can be directly observed from Fig.3. The versor differs from 
the two remaining by that it is parallel to BP and waterline — 
— therefore parallel to common edge of both planes, i.e the 


e = (cosy , siny , 0) 
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trace of BP on waterline. Let's observe that for œ = 0 the rota- 
tion axis e is undefined as the above mentioned planes are pa- 
rallel and their common edge does not exists; the product kxn 
is also undefined as both versors are parallel. 

As rotation axis in ship's upright position is undefined its 
choice is arbitrary. It is usually assumed that for & = 0 the an- 
gle y = 0, i.e. that in the initial position the rotation axis coin- 
cides with the trace of PS on waterline, however not necessari- 
ly. As the rotation axis an arbitrary axis on waterline can be 
chosen, especially its main inertia axis. From the previous con- 
siderations on the floatation axis it follows that in the latter 
case it must permanently coincide with the rotation axis (the 
angle x = 0, i.e. f= e) for every angle of heel Eq. (36) is valid 
regardless of the choice of rotation axis for œ = 0 (an axis fixed 
in space). 

The rotation axis is determined by a direction of heeling 
moment action in a ship-fixed reference system. As for the 
same analytical angles ọ and 9 the rotation axes are differently 
located relative to ship depending on an assumed reference axis, 
the rotation planes are also different, and hence the righting 
arms are different. This can be clearly stated from Eq.(10) where 
the righting arm GZ depends on r , e and n. For the same an- 
gles @ and 9, the versor n is the same but the rotation planes 
have different e and r, which makes GZ arms different. More- 
over if in physical space the rotation planes are the same as in 
Fig.2 , 7 and 8, rotated by the same angle n, then the analytical 
angles are different, which gives different curves of righting 
arms in function of the rotation angle n taken as ship's heel 
angle. 


ha 


e(x 


BẸ 
Fig. 9. Top view of the ship heeled by the angle ġ = 90° 


The case of the ship heeled by 90° is interesting as diffe- 
rences between various reference axes can be then distinctly 
observed. In the case of the reference axis Oy at the angle ọ = 
= 90° PS is horizontal, and the rotation plane crosses ship's 
gravity centre and buoyancy centre (Fig.9), which means that 
the ship is longitudinally balanced. As the plane is motionless 
in space the entire figure should be horizontally rotated around 
the point G (in this case to the left) by the angle of inclination 
relative to the axis z. The righting arm is negative and equal to 
the horizontal distance between the points G and B. 

For the axes x and z the rotation planes in Fig.9 cover each 
other — as they cross the point G parallel to the axis z. As the 
buoyancy centre is not located on the plane, the ship would so 
trim (in this case by bow) as to make the centre to be located 
on the rotation plane. 

In both cases the righting arm would be the same, equal to 
the horizontal distance between the gravity centre and transla- 
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ted buoyancy centre. However the heel angle n, in the sense of 
the rotation angle of the plane of rotation, would be different. 
For the axis x the angle n will be less than 90°, even by 
a few degrees, if trim at large heel angles is significant, which 
is possible especially in the case of a damaged ship; and for the 
axis z the angle n = a = 90°. 

The previous considerations dealing with the metacentric 
radii, floatation axes and curve of righting arms do not lose 
anything in value if to replace the differential do with the diffe- 
rential dn, and to substitute n for p . The directed angle fda, 
is inclined by the angle % to the appropriate rotation axis e 
defining the axis € in Fig.4. The axial component of the angle 
dn = dacosy, is the rotation angle of the rotation plane. The 
component da normal to the rotation axis is given by the for- 
mula analogous to Eq. (23) : 


- dO, , for the reference axis Ox 
(37) 
- dy-sino, for the reference axis Oz 


do = do, ‘siny = dn-tgx = 


From the above given equation 
two conclusions analogous to the previous ones, result : 


> the more deflected the floatation axis from the rotation axis 
the greater trim changes during heeling, and 

> when x = 0, i.e. when the floatation axis is parallel to the 
rotation axis, the ship trim does not change due to rotation 
of the waterline itself. As opposed to the reference axis Oy, 
there is no problem with determining the rotation axis at 
90° heel for the remaining reference axes. 


In the case of the reference axis Ox the rotation axis of the 
midship section around the normal id@ has the horizontal com- 
ponent dn = d@cos®@, which is axial one. The angle id@ has 
also the vertical component equal to - d@sinO, = - dntg0x. 

In the case of the reference axis Oz the vertical component 
of the rotation angle of BP around the normal kdw, equals 
dycosc. The vertical component makes the angle % changing 
by d%o as a result of the change of orientation of rotation axis 
(relative to the hull), caused by the rotation of the waterline 
alone. The angle dy is equal to the vertical component of the 
rotation. Hence one obtains : 


- dy-tg0, , for the reference axis Ox 
dxo = (38) 


- da-tgy/tga , for the reference axis Oz 


In the second part of Eq.(38) the identity appearing in Eq. 
(37) was utilized. The remaining comments concerning equi- 
-volume waterlines are still valid also for other reference axes. 
In considerations dealing with cross-curves of stability the an- 
gle þ should be obviously replaced with a relevant angle n. 
The quantity Ax. = Azge,, where e; is z-component of the 
relevant rotation axis e. The differential quotient amounts to : 

Al/Azg = eztgx. 


RIGHTING ARM CURVE 
OF MINIMUM STABILITY 


As previously mentioned, most heeling moments affecting 
the ship is parallel to PS, therefore a free - floating ship takes 
such position relative to its rotation plane as to make trace of 
water on PS normal to it. Therefore the question arises which 
position the ship takes when direction of the moment is not 
associated with its orientation relative to ship's hull. In other 
words, which reference axis should be then chosen. In order to 
unambigously answer this question it is necessary to know the 
mechanism of taking position by the ship relative to its rota- 
tion plane in the case of a free-heeling moment. 


In the case of a heeling moment not associated with the 
hull, such, for instance, as that resulting from shifting a weight 
on ship, or taking a weight to any place on ship, the ship in 
question rotates in such a way as to make potential energy of 
the heeled ship minimum [13]. In other words, the ship heels 
in such a way as the work associated with its heel is minimum. 
As moment work is proportional to the dynamic arm the mini- 
mum value of potential energy corresponds with minimum of 
the dynamic arm. From the classical ship theory it is known 
that dynamic arm depends on a course of metacentric radii in 
function of the heel angle, which — for a free-floating ship — 
— means : in function of the rotation angle of the rotation plane. 


Hence, in general case : 


la = J] BMsin m - v)dv -a (1 -cosn) (39) 
0 
where : 
v - dumb variable of integration, varying from 0 to n 


(the assumed rotation angle of rotation plane) 
BM - metacentric radius given by Eq.(12) 
a=BoG - constant value. 


Obviously the integrand in Eq. (39) decides on minimum 
dynamic arms, which takes its miniumum value for minimum 
metacentric radii in function of rotation angle. And these appe- 
ar when instantaneous floatation axis is permanently parallel 
to the rotation axis (f = e), i.e. when the angle % = 0. Hence 
a free-floating ship under influence of a free-heeling moment 
takes such position relative to its rotation plane as to make the 
instantaneous axis of floatation normal to the rotation plane. 


A few conclusions yield from that : 


“* From Eq.(15) it results that the rotation axis coincides — 
— then and only then — with the floatation axis when the 
deviation moment D = 0. The floatation axis is then the 
main inertia axis of the waterplane, for which the transverse 
inertia moment is minimum. When the ship heel is being 
changed a new floatation axis rotates relative to the previous 
one by the angle d,-= dx, which results only from the change 
of orientation of the main inertia axis of waterline since 
dXo = 0 due to zero value of the angle %, which can be 
deduced from Eq.(38). 


As the floatation axis is permanently parallel to the rotation 
axis Eq.(11) is reduced to the form : dn = da, = da. And, 
this in turn means that the rotation angle of rotation plane is 
equal to the angle contained between the actual waterplane 
and the initial one. Hence the curve of righting arms is 
a function of the angle a, and Oz - axis is the apropriate 
reference axis. 


\? 


% 


2e 
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As the angle x is equal to zero for any angle œ the correc- 
tion Al takes zero value, which can be seen from Eq. (29). 
Therefore righting arm curve can be found on the basis of 
cross-curves of stability by neglecting the correction as in 
the case of even-keel ship. 


e 


% 


For a given angle it is necessary to find such value of the 
angle of rotation around the versor no, W , as to obtain the 
centre of buoyancy located in the rotation plane containing 
the ship centre of gravity, and perpendicular to the longitu- 
dinal axis of inertia of waterline. 


2e 


% 


Righting arms are then as small as possible and such righ- 
ting arm curve is called the curve of minimum stability. Any 
other curve can have the same or greater arms. 


Worth mentioning that NAPA software as well as many 
other computer programs used in West countries for ship stabi- 


lity calculations have no possibility to calculate the righting 
arm curve of minimum stability as such mode of calculations 
is there unknown at all. 


CONCLUSIONS 


This paper presents the theoretical background for deter- 
mining the righting arm curve for a free-floating ship, under 
assumption that the righting moment has constant direction in 
space. It is also highlighted what is meant by the heel angle 
and how ship's longitudinal balance is understood. Results of 
the performed considerations can be summarized as follows : 


> Mode of calculation and choice of a reference axis ( Euler 
angles) considerably influence the results of calculation 
of the righting arm curve, hence they should be obligatori- 
ly defined. Such recommendations are not included in 
the current rules, which sometimes leads to great discre- 
pances between calculation results obtained from different 
sources. 


> Righting arm curves should be calculated — as a rule — for 
a free-floating ship balanced in an appropriate rotation pla- 
ne fixed in space and passing through the ship centre of 
gravity. The righting arm curve is then a function of the 
rotation angle of the rotation plane around the normal axis 
(rotation axis). 


> The rotation plane of an intact ship is perpendicular to tra- 
ce of water on PS, being a nodal line,due to the assumed 
character of heeling moments resulting mainly from hori- 
zontal, wind-generated forces. A rotation angle is the slope 
angle of y-axis relative to water-level, and a trim angle is 
the slope angle of trace of water on PS relative to BP. They 
are Euler angles associated with y-axis. 


> For a damaged ship which mainly heels under influence of 
vertical forces the rotation plane is perpendicular to the main 
inertia axis of waterline, being a nodal line. A rotation an- 
gle is the angle & contained between the initial waterline 
and water-level, and a trim angle is the angle contained 
between trace of water-level and trace of PS on the initial 
waterline. They are Euler angles associated with the ship- 
-fixed axis normal to the initial waterline. The calculations 
yield the righting arm curve of minimum stability. 


> The notion of cross-curves of stability is valid also for 
a free-floating ship heeled around the main inertia axis of 
instantaneous waterline. In other cases the calculated righ- 
ting arm curve should be corrected by introducing the cor- 
rection which accounts for an influence of change of trim, 
resulting from change of height of ship's centre of gravity. 


> It is advisable to perform calculations of the righting arm 
curve by means of the equi-volume waterplane method 
(Krilov-Dargnies's) which in natural way tracks translation 
of buoyancy centre during heeling [4]. The method has not 
been used so far in calculation practice, though it appears 
not only more accurate and well numerically conditioned, 
but also it significantly shortens time of calculations (about 
20 times) in comparison with bouyancy methods. Therefo- 
re it is worth implementing into common practice. 
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NOMENCLTURE 
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— height of gravity centre over buoyancy centre in upright 
position of ship 

— buoyancy centre 

— transverse metacentric radius 

— longitudinal metacentric radius 

— height of gravity centre over buoyancy centre 

— direction of rotation axis (unit vector normal to rotation 
plane) 

— versors (unit vectors) of trace of water on PS and on 
midship plane, respectively 

— versor (unit vector) of floatation axis 

— freeboard 

— gravity acceleration 

— ship gravity centre 

metacentric height 

righting arm 

versors (unit vectors) of the ship-fixed reference system 

whose origin is in the point K (intersection point of the 

plane of symmetry, PS, midship plane and base plane BP) 

— righting arm and dynamic arm, respectively 

length, breadth and mean draught of ship, respectively 

— upward pointing versor (unit vector) nomal to waterline 

— volumetric displacement of ship 

— angle between initial waterline and water-level 

— angle between trace of water on PS and midship plane 

— ship buoyancy (weight of displaced water) 

— correction of righting arm obtained by means of cross- 
curves of stability, accounting for oblique translation of 
gravity centre relative to rotation plane, due to changing 
the height of ship gravity centre over BP 

— rotation angle of plane of rotation 

— slope angle of trace of water on PS relative to x-axis 
of ship 

— water density 

— slope angle of trace of water on the stations relative to 
y-axis of ship 

— angle of PS inclination from the vertical 

— angle between floatation axis and rotation axis 

— angle between water-level trace line and PS trace line on 
initial waterline. 


€ 


Footnotes 


BoYZ reference system is fixed to the rotation plane and its ori- 
gin is located in the initial position of the buoyancy centre Bo. 
e.g. such as STATAW and SEA software systems used by CTO 
and PRS. 

In general case the nodal plane is a plane parallel to waterline of 
ship in the upright position, which rotates together with it. Its 
normal versor Ng is the same as that of initial waterline. The angle 
& is determined by : cos& = non, and y is the angle between 
trace of water and trace of PS in nodal plane. 


Acronyms 


BP -ship's base plane 

PS -ship's plane of symmetry 

CTO - Ship Design & Research Centre, Gdansk 
PRS -Polish Register of Shipping 
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A parametric method for preliminary 
determining of mass characteristics 
of inland navigation ships 


Jan P. Michalski 
Gdansk University of Technology 


ABSTRACT 


This paper presents a method for estimation of mass characteristics of vessels, elaborated 


with the use of an algorithm based on requirements of the Rules for the Classification and 
Construction of Inland Waterways Vessels of Polish Register of Shipping, and on a simpli- 
fied method [6] for determination of mass of hull plating stiffeners. The dimensioning 
method of hull structure scantlings based on this algorithm concerns classical vessels 
intended for the carrying of general cargo, dry and liquid bulk cargoes, which determine 
the range of the method application. The method does not cover vessels of different con- 
struction, e.g. roll-on-roll-off type vessels equipped with heavy decks, as the dimensioning 


of their scantlings is based on different relationships and models. 


Keywords : Inland waterways vessels, ship preliminary design methods 


INTRODUCTION 


Methods for preliminary, approximate estimation of mass 
as well as mass centre coordinates of a designed ship is an 
important branch of ship design theory. Knowledge of the pa- 
rameters is necessary already in early design phases — to itera- 
tively balance ship's floatability and stability before structural 
strength calculations are performed. 

The subject-matter literature contains some methods con- 
cerning this subject, e.g. in [1, 2, 3, 4] and [5], which, however, 
deal with sea-going ships whose features are different from 
those of inland waterways vessels. The methods applicable to 
inland waterways vessels are rather scarce, e.g. [6, 7, 8, 9, 10]. 

In the preliminary design stage, ship mass and its centre of 
gravity are usually estimated by means of the methods making 
use of parent ship's data, if only they are available. In the case 
when such data are lacking, or if a design project covers a bro- 
ad range of variability of ship design parameters, e.g. in a de- 
sign optimization study (with a selected objective function), 
then use of general parametric methods is necessary to deter- 
mine ships mass characteristics. 

The parametric methods for estimating ship hull mass are 
elaborated on the basis of either statistical data of existing ships 
or the rules of classification societies. 


A range of applicability of such methods 
is usually limited by the following factors : 


> ship’s functional type 

> hull structural arrangement 
> structural material 

> ship’s size range. 


To elaborate the presented method a set of discrete values 
of mass of vessel's hull series of systematically varying para- 
meters was simulated by means of the computer implementa- 
tion [14] of the algorithm [13]. Next, correlations between the 
hull parameters and their mass characteristics were investiga- 
ted, which made it possible to reduce some number of unimpor- 
tant parameters appearing in the problem. Next, analytical rela- 
tionships approximating the discrete values of mass characteri- 
stics were determined. Simulating calculations were carried out 
for assumed parameters of normal strength (NW) steel — 
— usually applied for construction of hulls of inland waterways 
vessels. 

The problem is to determine an analytical mathematical 
model to transform those vessel's features which are known in 
the preliminary design stage, i.e. 


@ the numerical parameters of ship x = (L, B, H, T, Cp) 
which describe : main dimensions, hull form, coefficients, 
structural material properties 

@ the qualitative atributes q = (K, R, Pi- Pa) which identi- 
fy vessel's functional type, class and topological features 
of hull structure 

@ the strength requirements of classification rules ; 


into estimated values of the hull mass M 
and its gravity centre height, Zg : 


M=F(x,q) Zg=F(X,q) 0 


A practical aspect of this work is to get a parametric me- 
thod easily codable into a computer program — useful for the 
estimating of the mass and gravity centre height of inland wa- 


and 
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terways vessels in the preliminary design stage. The cognitive 
aspect of the work is to investigate and broaden theoretical 
knowledge on the relationships between parameters of inland 
waterways vessels and their mass characteristics. 

The research in question has resulted from practical needs 
which have arisen during studies on the INCOWATRANS 
E!3065 (Inland and Costal Water Transport System) Project 
carried out within the frame of the UE EUREKA Program 
aimed at development of a modern, ecological friendly fleet of 
inland waterways vessels. 


RESEARCH PROGRAM 
AND ITS ASSUMPTIONS 


By making use of relevant knowledge dealing with sea- 
-going ships, given in [5,15] the hull parameters a priori assu- 
med to be significantly affecting its mass characteristics, were 
preliminarily selected. Hull mass magnitude depends on both — 
— vessel's numerical parameters and its qualitative features. 
Some assumptions simplifying real relationships were introdu- 
ced to determine simple mathematical models — describing in 
an approximate way the problem in question. It was assumed 
that for purposes of approximate estimation of vessel's mass its 
hull can be identified enough correctly by the following items : 


two functional types of vessel : 
dry cargo vessels 
liquid cargo vessels 


vessel’s class 
defined by permissible region of navigation 


vessel’s main particulars : 


L - length 

B - breadth 

H - depth to upper deck 
T - design draught 

Cg- hull block coefficient 


other parameters : 
k - permissible stress coefficient 
Re - yield strength of structural steel. 


The hull structure scantlings of dry cargo vessels and li- 
quid cargo ones are determined in a different way — therefore 
the following structural elements are taken into account [13] : 


for dry cargo vessels : 
¢ hatch openings in decks 
+ hatch coamings 
* specific features of double bottom 


for liquid cargo vessels : 
¢ longitudinal bulkheads 
* continuous structure of decks 
* specific features of double bottom. 


Variability ranges of the hull parameters of investigated hull 
series were determined on the basis of the data from literature 
[10, 17, 18]. Calculations for hull series covered all combina- 
tions of the following items : 


* functional types of liquid and dry cargo vessels 

* PRS (Polish Register of Shipping) ship classes concerning 
the restricted navigation regions R : 1, 2, 3 

* values of the length L : 25.00 m, 50.00 m, 75.00 m and 
100.00 m 

* values of the breadth B — resulting from the assumed ratio 
of L/B : 4, 5,6 

* values of the depth H resulting from the assumed ratio of 

L/H : 10, 15, 20 
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* values of the draught T resulting from the assumed ratio of 
H/T : 1.25, 1.5, 2.0. 


Moreover, the following constant parameters were assumed : 


% values of the block coefficient Cp : 0.7, 0.8, 0.9 
% normal strength (NW) steel of 235 MPa yield strength 
% permissible stress coefficient k = 0.7. 


The number of hull variants I,, of the hull series is equal to 
product of variants of values of varying parameters, i.e. : ves- 
sel class, length, breadth, depth, draught and block coefficient, 
namely Iw = 972 variants of hulls of any considered vessel 
functional type. 

Simulating calculations were performed by using the 
elaborated computer program INLAND VESSEL _ HULL_ 
_MASS.PAS [14] which realizes the algorithm for determi- 
ning hull scantlings of inland waterways vessels [13]. Results 
of the calculations are presented in the report [14]. 


MODELING VESSEL'S 
MASS CHARACTERISTICS 


An analytical relationship expressing the relation between 
mass of hull and its main design parameters should have the 
following features : 


* accuracy — to correctly approximate a statistical sample 

* simplicity — mathematical model should be convenient for 
design calculations and computer programming 

* good estimating properties — to provide credible predic- 
tions. 


It is not an easy task to satisfy the above mentioned requi- 
rements because of multi-dimensionality and nonlinear cha- 
racter of the analyzed relationships. To find the relations the 
least square approximation method was applied, and by means 
of the iterative procedures : 


+ hypotheses on approximation relationships were stated 

+ structural parameters of approximation model were deter- 
mined 

+ preliminary tests were performed 

+ approximation hypotheses were verified by assessing mea- 

sures of accuracy. 


To obtain sufficient accuracy a special normalizing techni- 
que was applied, similar to that given in [15]. The procedure 
consisted in preliminary choice of the normalizing function 
g(X) monotonically compatible with the approximated func- 
tion F(x) -— representing discrete values of simulated masses of 
hull series , and in determining the normalized smooth func- 
tion W(X) (of small variations) : 


F(x = 
© mit) 
g(x) 
Next, the function a(x) approximating the normalized 
function w(X) is determined. If the determined function a(x) 


correctly approximates the function w(x), then the function 
f(x): 


XEQ 


(2) 


f(x) = g(x)-a(x) 


correctly approximates the function F(x), hence this is the se- 
arched for solution of the problem in question. 

In the case of modelling hull mass characteristics a norma- 
lizing function may be e.g. either the product of vessel main 
dimensions (modular function) : 


e(x)=L-B-H 


(3) 


(4) 


or the power function : 
g(x) =c-] [ x# =c- L" B? H” -T -C R" (5) 
i=0 


The proper selection of a normalizing function is important 
in order to get sufficient accuracy of hull mass characteristics 
approximation. 


APPROXIMATED HULL MASS 
CHARACTERISTICS 


Final results of the best parametric models for mass cha- 
racteristics are below presented in the tables containing the 
following items : 


%* hypothesis on a form of normalizing function 

æ hypothesis on a form of approximating function with a list 
of variable parameters (great-letter symbols) and structural 
constants of the model (small-letter symbols) 

%* estimated values of structural constants of the mathemati- 
cal model 

#* correlation indices for assessing significance of design pa- 
rameters 

% estimated values of approximation accuracy measures, 
where : 


e E-average [%] — relative percentage error of approxi- 
mation 

e E-max [%]— maximum relative percentage error of ap- 
proximation for the whole sample 

e E >x% [%] — percentage of relative errors exceeding 
x% in the whole sample (of 972 elements). 


The requirement for model simplicity was fulfilled by ap- 
plying the power function formulas only — useful in calcula- 
ting values of the function derivatives — a feature convenient 
for linearization of a mathematical model of a design problem. 


Hull structure mass estimation 


Significance of the influence of vessel parameters on hull 
structure mass was assessed by means of the determined coef- 
ficients of correlation, then non-significant parameters — of 
small correlation coefficient values — were eliminated from the 
approximation model. In the case of liquid cargo vessels the 
obtained best approximation formula is given in Tab. 1 


Tab. 1. Structure, constants and accuracy 
of the approximation formula for hull mass estimation 


Liquid cargo inland vessels — sample of 972 simulated hulls 
Approximation hypothesis : g(x) =L-B-H 


M,,(&) — g(x): f(x) = L” : Be . H” : T* . Ce Ro 


Hypothesis verification 
Determined values of structural parameters 


| e | 0.0631615| et | 1.43625 | eb | 0.973853 | 
et | 0.150033 | ce | 0.154456 | er | -0.264760 | 


fe nee see |e 
E> 10% (a 


In the case of dry cargo vessels the best approximation was 
obtained by the formula whose description and accuracy as- 
sessment is given in Tab.2. 


Tab. 2. Structure, constants and accuracy 
of the approximation formula for hull mass estimation 


Dry cargo inland vessels — sample of 972 simulated hulls 
Approximation hypothesis : g(x) =L-B-H 


M,, (X) = g(x): (f, (&)+ £,(%)) =F, (&)+ F, (&) 
FE) =¢-L" > B”. H” T*<C,” R" 
LsL B HTAR" 


Hypothesis verification 
Determined values of structural parameters 


zoss [sp isan [oh | 14207 
C [sr [ose |_| 


Approximation accuracy measures 
E-average [%] E-max [%] E> 10% [%] 
10.08 


Estimation of height 
of hull structure gravity centre 


No correlation was revealed between the gravity centre 
height and the hull block coefficient Cp as well as the naviga- 
tion region R — both for dry cargo and liquid cargo vessels. 

The best approximation model found for the gravity cen- 
tre height of liquid cargo vessels, Zg =f (x), is presented in 
Tab.3 and in the case of dry cargo vessels — in Tab.4. 


Tab. 3. Structure, constants and accuracy 
of the approximation formula for hull gravity centre height 


Liquid cargo inland vessels — sample of 972 simulated hulls 
Approximation hypothesis : g(x) =H 


Zg(x) = g(x)-f(x)=c-L" -B?-H"-T" 


Hypothesis verification 


Determined values of structural parameters 
Pe [0.46023 [et] 0.00939] eb | 1. 7816-5 [ch [09250 [et] 275365 


Variables] L | B | n | T | œf R| 
Ze | 0826 | 076s | 1000 | 0932 | o.000 | 0.000 
E-average [%] E> 5% [%] 

Do oeo |o s | oo 


Tab. 4. Structure, constants and accuracy 
of the approximation formula for hull gravity centre height 


Dry cargo inland vessels — sample of 972 simulated hulls 
Approximation hypothesis : g(x) =H 


Z9(X) = g(x)-f()=c-L" -B®-H" -T* 


Hypothesis verification 
Determined values of structural parameters 
| e | 0.40144 | et | 0.09297 | eb | -1.960E-5 0.7177 | et |0.000307 


Variables] L | B | n | T | G | R | 


E-average [%] 
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E-max [%] 


E > 5% [%] 
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POWER PLANT MASS 
AND OUTFIT MASS CHARACTERISTICS 


Estimation of mass of medium 
and high-speed marine diesel engines 


Ship power plant mass essentially depends on a type 
of engines of the propulsion system. On inland waterways ves- 
sels, are commonly installed medium-speed diesel engines 
(of 500+600 rpm) and high-speed ones (of 1500+1800 rpm), 
which transfer power either directly, or through mechanical, 
hydraulic or electrical transmission gear — to the screw pro- 
peller(s) of fixed or controllable pitch, often accommodated 
in a Kort nozzle. 

Diesel engine mass depends on its rated torque, power out- 
put, cylinder arrangement and its supercharging. Apart from 
the engine itself, mass of power plant consists of masses of its 
equipment and systems. Data contained in the subject-matter 
literature usually concern ship power plants with the engines 
of low speed — applicable to sea-going ships; e.g. in [1] is gi- 
ven a formula for determining mass of ship diesel engines of 
low and medium speed — of more than 1000 kW rated power : 


P 0.84 
Min = fz | [t] (6) 


n 
where : 
P - installed engine power [kW] 
n - engine speed [rpm]. 

Publications dealing with vessel power plants equipped with 
medium-speed and high-speed diesel engines of relatively small 
power are scarce. 

For purposes of the INCOWATRANS design study on in- 
land waterways passenger vessel, a research work on mass of 
medium-speed and high-speed engines of small power was 
performed on the basis of the catalogue data of: Volvo-Penta, 
Deutz, MAN, Caterpillar, Detroit Diesel, Wola-Henschel and 
Wartsila. The collected data of the engine parameters and its 
mass values were used as a statistical sample to determine an 
analytical relationship approximating mass of medium-speed 
and high-speed marine diesel engines of small power. In re- 


sult, the following simple formula was obtained : 
0.85 


Min STIN gs. [t] (7) 


where : 
P - installed engine power [kW] 
n - engine speed [rpm] 
A comparison of the unit mass of the engines : 


Muy [kg/kW] 


Me = 


calculated by the formula (6) — dashed lines (engines for sea- 
-going ships), and (7) — continuous lines (engines for inland 
waterways vessels) is presented on Fig.1. 


From the comparison it results that in the case of medium- 
-speed engines (of 750 rpm ) the engine mass according to (7) 
is greater by about 60% than that estimated by using (6). For 
the high-speed engines (of 2500 rpm) the mass according to (7) 
is smaller by about 50% that that estimated by means of (6). 
For the engines of speeds within the range of 1000+2000 rpm 
so estimated values of mass are similar. 


Estimation of mass of power plants fitted with 
medium and high-speed engines 


Power plant mass of sea-going ships without stores (,,dry 
weight” condition) acc. [1] can be determined as a sum of mass 
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Watson [1] 


Unit mass m, [kg/kW] 


0 400 800 1200 


Installed engine power P [kW] 


1800 2000 


Fig. 1. Comparison of the values of the unit mass of ship medium- and 
high-speed diesel engines, obtained from the formula (6) (dashed lines) 
and (7) (continuous lines) 


of the main engine and mass of the remaining elements of the 
power plant : 


M, =k: P” [t (8) 


where : k, - a coefficient taking values from the interval 
of 0.19 + 0.83 - depending on a ship functional type. 


Main engine mass of inland waterways vessels reaches 
about 30% + 45% of the total mass of the power plant (in "dry 
weight" condition). To determine mass of a vessel power plant 
fitted with medium-speed or high-speed engines the following 
formula was elaborated : 

0.85 


M, =7130:— +k,:P””’ [t] (9) 
zi 


Fig.2. presents a comparison of the unit mass of power plants 
Mp 
m, = P [kg/kW] 


calculated by means of the formula (6) and (8) — dashed lines, 
and formula (9) — continuous lines — in both cases kg = 0.2 


[kg/kW] 


acc. to (9) 


P 


acc. to (6) and (8) 


Unit mass m 


800 1200 
Installed engine power P [kW] 


1800 


Fig. 2. Comparison of the unit mass of ship power plants with medium- and 
high-speed diesel engines, obtained from the formula (6 and 8) — dashed 
lines, and formula (9) — continuous line 


In the case of power plant of high-speed engines it is often 
necessary to apply a reduction gear. To estimate mass of 


a mechanical reduction gear the following formula was elabo- 
rated : P 
M,=4-— i (10) 
n 


Estimation of hull outfit mass 


Mass characteristics of hull outfit are not easily modeled 
mathematically as they depend not only on numerical parame- 
ters of a ship but also on such factors as: ship functional type, 
outfitting standard, crew number, as well as on obligatory re- 
quirements of relevant provisions of maritime administration 
and a society classification rules. Within the frame of the IN- 
COWATRANS project no such research work was made. But 
on the basis of the diagrams given in [7] the following appro- 
ximating analytical formula concerning inland waterways ves- 
sels, useful in the preliminary design stage, were elaborated : 


hull equipment mass : 


M, =0.125-(L-B-H)*” [t] (11) 
mass of hull systems : 
M, =0.25-(L-B-H)'” [t (12) 
hull outfit mass : 
M, =0.025-(L-B-H)** [t (13) 


where all dimensions are in [m]. 


PRELIMINARY ASSESSMENT 
OF THE METHOD'S 
PREDICTING CAPABILITY 


An attempt to preliminarily assess predicting capability of 
the method in question was made by taking into account, as an 
example, several vessels whose data were obtained from their 
PRS class documentation. In the case of the SINE 207 inland 
waterways container vessel of 18/36 TEU capacity its prelimi- 
nary design documentation was used. Results of the verifying 
calculations are presented in Tab.5. 


SUMMARY 


O The scantling determination algorithm used to elaborate the 
presented parametric method concerns the inland waterways 
vessels intended for the carrying of general cargo and dry 
and liquid bulk cargoes. 


O The performed verification of the method showed that in 
the case of the passenger vessels and container vessels the 
achieved predictions of : 

* mass characteristics — are acceptable approximations 
for purposes of the preliminary design stage (mean er- 
ror is of about 15%) 

* gravity centre height characteristics — are acceptable 
approximations for purpopses of the preliminary design 
stage (mean error is less than 14%) 


O In the case of the vessels of a different structure, e.g. car 
ferries — intended for the carrying roll-on — roll-off cargoes, 
fitted with heavy decks (M/V Berlin, M/V Bielik ID), the 
obtained results are different from the values summed up 
of detail mass specifications. 


Q For light cargo vessels the proposed parametric method of mass 
prediction can be a useful tool in preliminary design stage. 


NOMENCLATURE 

B - breadth 

Cg - block coefficient 

E - relative percentage error 

f(X) - approximating function of hull mass and its gravity centre 

F(X) - approximated function of hull mass and its gravity centre 
height 

g(x) - normalizing function 

H - depth 

L - length 

m, - unit mass of engine 

m, - unit mass of power plant 

Me - hull equipment mass 

My - hull mass 

M, - real vessel mass 


Tab. 5. Results of verifying calculations to assess predicting capability of the presented method 


Parameters 


Symbol 
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NAVAL ARCHITECTURE 


NAVAL ARCHITECTURE 


Mrp - predicted mass of vessel 

Mo - hull outfit mass 

Mp - power plant mass 

Ms - mass of hull systems 

n  - engine speed 

NW - normal strength steel acc. PRS 

P - installed (rated) engine power 

q  - vector of qualitative design attributes of vessel 
R - permissible region of navigation 

T  - draught 

x - vector of numerical design parameters of vessel 
Z, - predicted height of vessel's gravity centre 

Zor - real height of vessel's gravity centre 

Q - set of acceptable parameters. 
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Ship Stability Workshop 


After the preceding workshops of the kind held in Great 
Britain, Japan, Greece, Canada, Italy and USA, the 7th 
International Ship Stability Workshop was organized in 
China. The interesting scientific meeting of naval archi- 
tects was held on 1+3 November 2004 in Shanghai Jiao 
Tong University. 


Its program contained 28 papers 
presented during the following sessions : 


Theoretical Development in Damage Stability 
Damage Survivability Assessment 
Assessment of Ship Stability Safety 
Theoretical Prediction of Intact Stability 
Experimental Investigation of Intact Stability 
Ship Dynamics with Water on Deck 

and Extreme Waves 

Stability Research in China 


oO 008000080 


Authors of the papers represented scientific centres of 
Canada, China, Finland, Germany, Greece, Holland, Ita- 
ly, Japan, Sweden, United Kingdom and USA. 

Among them it was Prof. M.Pawtowski from Gdansk 
University of Technology, Poland, who took part in the 
activity of International Standing Committee, and — du- 
ring 2nd session — presented the paper titled : 


General framework of new subdivision regulations 


ce : Beata Gis 


Miscellanea 


Przemystu Okretowego 
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A new handbook written by Prof. B. K. Mazurkiewicz 
has been recently published (in Polish) by the Foundation 
for the Promotion of Maritime Industry, Gdansk. 


The book consists of five following chapters : 


1. General information 
This chapter (of 46 pages) contains general information 
on yachts, definition and classification of yacht ports as well 
as information dealing with their development in the light 
of the integrated management of coastal areas. Moreover 
prospects for building and modernizing yacht ports in Po- 
land are also described. 


2. Location, general arrangement 

In this chapter (of 51 pages) , after presentation of condi- 
tions for location of yacht ports, with special attention paid 
to those demanded by the environment, are presented infor- 
mation and recommendations concerning elaboration of ge- 
neral arrangement and spatial plan of a yacht port - marina. 
Moreover, parameters of the yacht port - marina as well as of 
inland waterways and water areas are described in detail. 


3. Solutions of hydro-engineering structures 
of yacht ports - marinas 
This chapter, the largest (of 136 pages) is devoted to the 
designing and construction of different hydro-engineering 
harbour structures such as breakwaters, quays, fixed and flo- 
ating platforms, mooring posts as well as locks and slipways. 
Especially great attention is payed to the breakwaters and 


Design 
of yacht ports - marinas 


for the Promotion 
of Maritime Industry 


==S——_—— 
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floating platforms, and first of all to anchoring such struc- 
tures. There are also given calculation methods mainly for 
floating platforms more and more frequently used in yacht 
ports - marinas. Due to the presented broad information on 
different floating platforms it is possible to make direct use 
of it in selecting parameters and solutions of the platforms 
functioning both as mooring ones and pedestrian passages. 


4. Equipment of yacht ports - marinas 

This chapter (of 45 pages) contains information and re- 
commendations concerning devices of different types re- 
quired for a yacht port - marina. These are : mooring devi- 
ces and fenders, electrical installations, water supply, fire- 
-fighting and telecommunication systems, devices and in- 
stallations for receiving, collecting and carrying away wa- 
ste and sewage, fuel supply stations for yachts, facilities for 
hoisting, transporting and launching floating units being re- 
paired or stored (wintering) in a given port as well as faci- 
lities making it possible to properly manage a yacht port - 
- marina and to ensure safety to floating units lying in it. To 
it belongs also navigation marking which makes it possible 
to safely reach an accomodation or repair berth. Summing 
up, the chapter covers the main requirements for design and 
use of different facilities and equipment necessary in a mo- 
dern yacht port - marina. 


5. Design of yacht port - marina 

In the 5th chapter (of 15 pages) is presented the scope 
of design of yacht port - marina with drawing attention to 
an appropriate sequence of design tasks to be elaborated, 
with reference given to particular chapters of the book. It is 
assumed that elaboration of design of a yacht port - marina 
on the basis of the proposed scope and course of designing 
would allow to avoid mistakes which usually demonstrate 
as late as during operation of a yacht port - marina. The 
chapter is illustrated by selected designs of yacht ports - 
- marinas. 


At the end of the book the used nomenclature and 
referred to bibliography (containing 81 items) is attached. 


iT 


Yacht port - marina in Leba (photo: C. Spigarski) 
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Experimental values of temperature distribution 
in a Sliding bearing sleeve 
lubricated with non-Newtonian oils 


Andrzej Miszczak 
Gdynia Maritime University 
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ABSTRACT 


This paper presents results of experimental investigations of temperature distributions on 
inner surface of a sleeve of transverse sliding bearing lubricated with non-Newtonian oils. 
The measurements were performed by means of Pt100 miniature sensors placed close to 
internal surface of the sleeve. To measure sleeve temperature distributions use was made 
of a test stand installed at Gdynia Maritime University. Delo” 1000 Marine 30 oil, SAE 
15W40 basic oil, and a ferro-oil made of the SAE 15W40 basic oil were used as a lubrica- 
ting medium. During measuring temperature distributions for the ferro-oil different 
intensity values of external magnetic field generated by electromagnets were applied. 


Key words : ferro-oil, non-Newtonian oil, temperature distribution, test stand 


AIM OF MEASUREMENTS 
AND DESCRIPTION OF TEST STAND 


The main aim of the measuring of temperature distribu- 
tion on the inner surface of the sleeve was to investigate influ- 
ence of magnetic particles contained in oil as well as of exter- 
nal magnetic field intensity values on temperature distribution 
in slide bearing. 

The experimental investigations was performed on the test 
stand installed at Gdynia Maritime University. The stand was 
built with making use of a part of T - 05 tester elaborated by 
Institute of Operation Technology, Radom (Fig.1.) The stand’s 
design was made on the basis of own observations of this au- 
thor and other designs presented in [1,2]. In the stand a new 
shaft with 100 mm diameter of journal was applied. The exis- 
ting driving motor was replaced with a new one of 2.2 kW 
output and twice greater rotational speed than original.The 
SJ100-022 HFE Hitachi inverter was applied to made continu- 
ous control of rotational speed possible. 

The stand in question makes also possible to measure rela- 
tive eccentricity and friction force at given operational para- 
meters such as : inlet oil pressure, load, rotational speed, tem- 
perature of sleeve or oil at inlet [3,4]. 

To measure friction force a tensometric gauge fitted with 
KT 1400/K/200N/2410D transducer of Megatron, was applied. 
Its measuring range for friction force reaches 200N at 0.17% 
non-linearity and 0.1N resolution, which gives 5mV output 
voltage. The gauge is supplied with 24V constant voltage and 
an output voltage signal is contained within 0+10V range, equi- 
valent to measured forces within 0+200N range. 

To determine relative eccentricity, two MDCT2/K/2410 
induction extensometers of Megatron were added to the stand. 
The extensometers had built-in transducers which made it pos- 
sible to read output signal in the form of standard voltage sig- 
nal of 0=10V range. 
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One of the extensometers was installed horizontally, the 
other — vertically [3,4]. They were of +1mm range at lum re- 
solution, which is equivalent to 5mV output voltage. When the 
measurements of the maximum sleeve translation to the left 
and right are made by means of the horizontal extensometer, 
and next the measurements of the maximum sleeve translation 
up and down with the use of the vertical one, then from the 
obtained results a value of the radial clearence of the sleeve, €, 
location of the journal centre O and that of the sleeve O’ , can 
be determined both in the hot and cold condition of the test 
stand. Additionally MDS10P Technicad contactless eddy-cur- 
rent sensor of journal’s translation was installed to take into 
account journal’s deflection under load. 


translation 
sensors 


Gy 


Friction force 
sensor 
L Eddy-current journal’s 
translation sensor 


Fig.1. General view of the test stand for determining temperature 
distributions, measuring friction force and relative eccentricity 
in journal sliding bearing 
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An additional oil tank fitted with a water cooler and 1 kW 
electric heater fed from an autotransformator was applied to 
deliver lubricating oil of appropriate temperature values to the 
sliding friction node. The delivered oil temperature was so con- 
trolled as to maintain it constant at oil inlet. In the lubricating 
system a PZ3-4/16-1-122 gear pump of 5+6 I/min capacity and 


6 bar maximum pressure, as well as a manometer and a set of 
control valves, were installed. 

The recording of indications of the translation transducers , 
of friction force transducer and temperature gauge , was per- 
formed by means of Catman computer software and UPM100 
100-channel signal recorder of Hottinger Baldwin Messtech- 
nik. With the use of the recorder standard voltage and current 
signals can be saved, temperature values can be measured by 
means of thermocouples and resistance gauges, e.g. of Pt100 
type, as well as indications of tensometers can be measured. 

A DC-supplied electromagnet (Fig.2) was built to investi- 
gate magnetic field influence on operational parameters and tem- 
perature distribution within the gap of the sliding bearing lubri- 
cated with ferro-oil. The maximum value of the magnetic field 
generated by the magnet was 100 mT. The magnetic saturation 
of the waste ferro-oil was 50+70 mT. Magnetic field values 
were measured with the use of a SMS-102 magnetic field me- 
ter fitted with a transverse sound, produced by Asonik. 


Mee 
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Fig. 2. General view of the applied electromagnet 


BEARING MATERIALS 


To determine temperature distribution over sleeve’s inner 
surface the sliding sleeve made of PBM 100120120 bronze by 
SKF, having the following dimensions : 100 mm inner diame- 
ter, 120 mm outer diameter and 120 mm length, was selected. 
In the sleeve several holes to accommodate Pt100 gauges in 
accordance with Fig.3 were made. The gauges were placed 


A 24° 


0.5 mm deep above the inner surface of the sliding sleeve. 
A way of measuring the temperature distribution was assumed 
the same as in the work [5]. Contact between a gauge and sle- 
eve was improved by covering the gauge with a heat-conduc- 
ting paste. The gauges were poured with an epoxy resin and 
connected with the UPM100 recorder by means of conductors. 
Fig.4 presents the arrangement of the gauges in the sleeve. 


Fig. 4. General view of the sliding sleeve 
for measuring temperature distribution and heat flux density 


Temperature distribution measurements were carried out 
for strictly determined operational parameters of the bearing 
and the following grades of engine oil applied as a lubricating 
medium : 


> Delo® 1000 Marine 30 TEXACO ship engine oil — in fresh 
condition i 

> the same Delo® 1000 Marine 30 — in a waste condition, i.e. 
after 300 h continuous operation within a ship engine (of 
an oil-rig support & salvage tug); the oil change time re- 
commended by ship’s owner was 300 h 

> Ferro-oil — made by FerroTec on the basis of SAE 15W40 
mineral oil free of any bettering additives (i.e. basic oil). In 
the oil were contained Fe;0, magnetic particles of 5+15 
nm in size and amount of abt. 1.5% vol., as well as a surfa- 
ce-active substance to prevent magnetic particles against 
clustering 

> SAE 15W40 basic oil — delivered by Gdansk Oil Refinery 
for testing 


Developed view of sliding sleeve 


Fig. 3. Schematic diagram of arrangement of Pt100 temperature gauges 
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TEMPERATURE DISTRIBUTION 
MEASUREMENTS 


Due to different properties of the tested oils as well as to be 
able to determine differences in results of measured tempera- 
tures, all the tested oils were subjected to the same selected 
operational parameters of the bearing. It was : its rotational 
speed, load, inlet pressure and temperature of lubricating oil. 
Tab.1. presents values of the operational parameters applied 
during measurements of temperature distribution. 


Tab.1. Values of operational parameters of tested bearing 


Inlet oil pressure [bar] 


Rotational speed [rpm] 


Load 
[N] 


D 


Inlet oil temperature : 70°C 


The measured radial clearance £ of the tested bearing was 
0.081 mm, which corresponded with the relative 
radial clearance y = 0.0016. 


0° and 360° 


a) eC 


Lubricant supply 


Lubricant supply 


120° 
240° 
360° @ 


Sleeve 


Lubricant supply 


Before the measurements the test stand was thermally 
stabilized by keeping a high speed of the journal for abt. 1h at 
simultaneous loading the sleeve with a weight. After such ope- 
ration performed in each case, appropriate values of the opera- 
tional parameters were set, some time was waited till stabiliza- 
tion of temperature indications takes place, and next the rele- 
vant temperature distribution measurements carried out. 

In this paper only some results selected out of a broad tes- 
ting program are presented in the form of diagrams of measu- 
red temperature distributions. 


Results for Delo® 1000 Marine 30 oil 
in fresh condition 


In Fig.5. are presented measured values of temperature dis- 
tribution over the surface of PBM 100120120 sleeve for Delo” 
1000 Marine 30 oil in fresh condition, under 1034.7N load, at 
1.5 bar pressure and its 70°C temperature of oil at inlet, as well 
as for three following rotational speeds : 840 rpm (Fig.5a); 
1960 rpm (Fig.5b); 2800 rpm (Fig.5c). In every figure are pre- 
sented the same two temperature distributions, but shown in 
a different axonometric view. 

9 


b/2 0° 


Fig. 5. Distributions of the temperature measured on the sleeve under 1034.7 N load at 1.5 bar pressure and 70°C temperature of oil at inlet, 
and the following values of journals rotational speed : a) 840 rpm; b) 1960 rpm; c) 2800 rpm 
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In all the diagrams, @ axes describe the direction of bearing rotational speed of the journal but two values of inlet oil pres- 
angle, beginning from the point of some degrees behind the sure. 
point of lubricating oil delivery, towards shaft’s rotation. The Fig.6a shows the temperature distribution measured on 
second horizontal axis marked z points along bearing’s length the inner surface of the sleeve at | bar inlet oil pressure, un- 
(b = 120 mm). In the middle length of the bearing a hole is der 1034.7 N load, and 70°C inlet oil temperature, and 2800 
made through which lubricating oil is delivered. The sleeve rpm rotational speed of the journal. Fig.6b shows the tempe- 
was not fitted with distinct oil pockets, but only with a small rature distribution for the same values of operational para- 
circular cavity abt. 0.8 mm deep, and of 10 mm diameter. meters of the bearing as for the Fig.6a, but for the different 

In Fig.6. is presented the influence of inlet pressure valu- value of inlet oil pressure, equal to 2 bar. The temperature 
es on temperature distribution in the same oil at only one distributions for 1.5 bar oil pressure are shown in Fig.Sc. 


Load = 1034.7 N 
Rotational speed = 2800 rpm 


Inlet oil pressure 1.0 bar 


9 


Sleeve 


0° and 360° [N 


360° © Lubricant supply 


Inlet oil pressure 2.0 bar 


Lubricant supply 


b/2 0° 


Fig. 6. Distributions of the temperature measured on the sleeve under 1034.7 N load 
and 2800 rpm rotational speed of the journal, for 70°C temperature of oil at inlet, 
and the different values of the oil pressure 
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Results for SAE 15W40 basic oil 


During tests for the SAE 15W40 basic oil the sleeve was subjected to the same operational conditions as those for the tests 
of the preceding oils. The relevant results of temperature distributions are shown in Fig.7 and 8. 


Load = 1034.7N Inlet oil pressure = 1.5 bar 


Rotational speed 840 rpm 
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0° and 360° (2 


Sleeve 
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Lubricant supply 
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Rotational speed 1960 rpm 


°C 


Lubricant supply 120° 
b/2 0° 


oÇ Rotational speed 2800 rpm 


240° 


Lubricant supply o 
360° Ọ 120 
b/2 0° 
Fig. 7. Distributions of the temperature measured on the sleeve under 1034.7 N load at 1.5 bar pressure and 70°C temperature of oil at inlet, 
and the different values of journal's rotational speed 
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Load = 1034.7 N 
Rotational speed = 2800 rpm 


Inlet oil pressure 1.0 bar 


0° and 360° 
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Fig. 8. Distributions of the temperature measured on the sleeve under 1034.7 N load 
and 2800 rpm rotational speed of the journal, for 70°C temperature of oil at inlet, 
and the different values of the inlet oil pressure 


Results for ferro-oil 


The measured temperature distributions concerned possi- 
ble changes resulting from changes of journal’s rotational spe- 
ed and simultaneous action of magnetic field induction at 
a given load exerted on the bearing. Three first figures present 
the temperature distribution for rotational speeds of : 840 rpm; 
1960 rpm; and 2800 rpm under 1034.7 N load on the sleeve 
and no external magnetic induction field. 


The measured values presented in Fig.10 were obtained for 
the same speeds, inlet oil pressure and load, but for the magne- 
tic induction field of about 30 mT intensity. 

Fig.11 presents the temperature distributions also for the 
same rotational speeds, 1034.7 N load and 1.5 bar inlet oil pres- 
sure, but for the magnetic induction field of about 55 mT in- 
tensity. 
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Ferro-oil without magnetic field 
Load = 1034.7 N Inlet oil pressure = 1.5 bar 


Rotational speed 840 rpm 
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Fig. 9. Distributions of the temperature measured on the sleeve under 1034.7 N load at 1.5 bar pressure and 70°C temperature of oil at inlet, 
and the different values of journal's rotational speed (no magnetic induction field) 
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Ferro-oil with magnetic field = 30 mT = 
Load = 1034.7 N Inlet oil pressure = 1.5 bar z 
oD 
Rotational speed 840 rpm a 
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Fig. 10. Distributions of the temperature measured on the sleeve under 1034.7N load at 1.5 bar pressure and 70°C temperature of oil at inlet, 
and the different values of journal's rotational speed (under magnetic induction field of 30 mT) 
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Ferro-oil with magnetic field = 55 mT 
Load = 1034.7N Inlet oil pressure = 1.5 bar 
Rotational speed 840 rpm 


0° and 360° 


Lubricant supply 
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Fig. 11. Distributions of the temperature measured on the sleeve under 1034.7N load at 1.5 bar pressure and 70°C temperature of oil at inlet, 
and the different values of journal's rotational speed (under magnetic induction field of 55 mT) 
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DISCUSSION OF THE TEST RESULTS 


Apart from the above presented results of temperature me- 
asurements also the Delo®1000 Marine 30 waste oil was 
investigated. The obtained results differed only slightly from 
those for the fresh oil of the same grade. Because of publi- 
shing limitations, results of temperature measurements for 
the remaining oils and other operational parameters, (e.g. 
for 2020 N load) are not included in this paper, neverthe- 
less they have been analyzed and discussed. 


Analyzing the values of temperature distribution one can 
state that the distributions for the fresh oil and waste one 
for the same operational parameters of the bearing differ 
only slightly. 


The rotational speed influence on the change of temperatu- 
re distribution in the Delo®1000 Marine 30 fresh oil, at 
1034.7N load, amounts to about 26% for the increase of 
journal’s rotational speed from 840 rpm to 1960 rpm, and 
for the increase from 840 rpm to 2800 rpm — to about 39%. 


In the case of the waste oil at 1034.7N load on the bearing 
the rotational speed changes influence on the change of 
temperature distribution amounts to about 24% for the in- 
crease of journal’s rotational speed from 840 rpm to 1960 
rpm, and to about 41% for the increase from 840 rpm to 
2800 rpm. 


In the case when the bearing is loaded with 2020N load, 
the influence of rotational speed on change of temperature 
distribution in the fresh Delo” 1000 Marine 30 oil amounts 
to about 23% for the increase of journal’s rotational speed 
from 840 rpm to 1960 rpm and to about 35% for the incre- 
ase from 840 rpm to 2800 rpm. 


In the case of the Delo®1000 Marine 30 waste oil and for 

2020N load exerted on the bearing the rotational speed in- 
fluence on change of temperature distribution amounts to 
about 24% for the increase of journal’s rotational speed from 
840 rpm to 1960 rpm, and to about 37% for the increase 
from 840 rpm to 2800 rpm. 


The influence of inlet oil pressure on temperature distribu- 
tion measured in the same points, both for the fresh and 
waste Delo” 1000 Marine 30 oil at 1034.7N load on the 
sleeve, is slight. The inlet oil pressure increased from 1 bar 
to 1.5 bar lowers temperature values by about 3%, and by 
about 6% in the case when the oil pressure is changed from 
1 bar to 2 bar. In the case when the sleeve is under 2020N 
load, the influence of inlet oil pressure on temperature dis- 
tribution is negligibly smaller than that of inlet oil pressure 
on temperature distribution at a lower load. 


Analyzing the results of measured temperature distributions 
in the SAE 15W40 basic oil one can draw similar conclu- 
sions as those for the above discussed oil, however with 
the difference that the temperature values for all the just 
presented cases are by several percent lower than those 
obtained for the previously tested oil. 


Hence the influence of rotational speed on change of tem- 

perature distribution in the basic oil at 1034.7 N load amo- 
unts to about 23% for the increase of journal’s rotational 
speed from 840 rpm to 1960 rpm, and to about 36% for the 
speed increase from 840 rpm to 2800 rpm. 


In the case when the bearing is under 2020 N load the rota- 
tional speed influence on change of temperature distribu- 
tion in the SAE 15W40 basic oil amounts to about 20% for 
the increase of journal’s speed from 840 rpm to 1960 rpm 


and about 30% for the speed increase from 840 rpm to 
2800 rpm. 


The influence of inlet oil pressure on temperature distribu- 
tion in the basic oil is also low. For 1034.7N load applied 
on the sleeve, and the same measurement points, it amo- 
unts to about 3.5% for the inlet oil pressure increased from 
1 bar to 1.5 bar, and to 6.5% — when the pressure is chan- 
ged from 1 bar to 2 bar. 


The influence of inlet oil pressure on temperature distribu- 
tion in SAE 15W40 basic oil under 2020 N load exerted on 
the sleeve, measured in the same measurement points, is 
equal to about 2% for the inlet oil pressure increased from 
1 bar to 1.5 bar, and 2.6% for the oil pressure increased 
from 1 bar to 2 bar, respectively. 


Analyzing the measured temperature distributions for fer- 
ro-oil, presented in Fig. 9+11, one can draw similar conclu- 
sions as those for the previously considered oils, with the 
difference that the temperature values for all the cases actu- 
ally considered are by a dozen or so percent greater than the 
temperature values obtained for the previously tested oils. 


When comparing temperatures of other oils with basic oil 
temperature it can be observed that the difference is even 
as great as 35%. Reminding that the ferro-oil was made of 
the basic oil, one can attribute its increased viscosity to the 
introduction of magnetic particles and the application of 
magnetic induction field. 


The influence of journal’s rotational speed changes on tem- 
perature distribution changes during lubricating the bearing 
with the ferro-oil was tested for 1034.7N load exerted on the 
sleeve. When the journal’s rotational speed was increased from 
840 rpm to 1960 rpm the temperature increased by about 29% 
(at OmT magnetic field intensity during measurements), 29% 
(at 30mT), 27% (at 55mT), and for the speed increase from 
840 rpm to 2800 rpm the temperature increased by 43% (at 
OmT), 42% (at 30mT), 40% (at 55mT) respectively. 


In the case when 2020N load was applied to the bearing the 
rotational speed influence on change of temperature distri- 
bution in the ferro-oil amounts to about 29% (at OmT); 24% 
(at 30mT); 23% (at 55mT) as a result of the increase of 
journal’s rotational speed from 840 rpm to 1960 rpm, and 
about 42% (at OmT); 36% (at 30mT); 35% (at 55mT) as 
a result of the increase of journal’s rotational speed from 
840 rpm to 2800 rpm. 


The influence of inlet oil pressure on measured temperatu- 
re distribution is, like for other oils, low and contained 
within the range of 2%+3% as a result of the change of 
inlet oil pressure within the range of 1+2 bar. 


CONCLUSIONS 


Considering the obtained test results both presented and 
not presented in this paper (due to a huge amount of them) 
one can draw the following conclusions : 


> The difference between the values of temperature in the 
tested waste oil, distributed on the inner surface of the 
sliding bearing sleeve in question and those in the fresh 
oil of the same grade in the same lubricating conditions, 
is small. The temperature increase by a few percent for 
the waste oil can be justified by its increased viscosity. 
The greater viscosity makes values of friction force in- 
creasing and thus more heat is emitted. At the same 
amount of the flowing oil and a greater amount of heat 
generated within oil film its temperature must be hi- 
gher than that of an oil of lower viscosity. 
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crease makes values of the sleeve’s temperature incre- 
asing , as well as the difference between the highest and 
lowest value of sleeve’s temperature greater. The high é 

temperature drops appearing in the central part of the I I 

diagram result from the delivering of lubricating oil at ls ce anea 
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> The distinct temperature drop resulting from the incre- 
ased load on the sliding bearing can be observed. The 
phenomenon could result from that at a high load the 
maximum lubricating gap increases, thus a greater amo- 
unt of lubricating oil could flow through the gap at an 
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A new method for searching optimal path 
on a raster plane including cost 
of direction changes 


Rafat Sztapczynski 
Gdansk University of Technology 


ABSTRACT 


The article introduces a new algorithm for finding optimal routes on raster planes. This 

method takes advantage of the new data structure and results in minimizing the number of 

direction changes within a route. It has linear time and space complexities and is suffi- 

ciently fast to perform real-time routing on the raster grids. Both the algorithm and its 

A data structure are presented in detail in the paper. Possible applications of this method 
i are also discussed. 


Keywords : navigation, optimal route, turn penalties, Lee’s algorithm, vave propagation, raster charts 


1. INTRODUCTION 


Raster grids are a digital representation of planary or spa- 
tial data, that is currently in use in many fields some of which 
belong to the maritime industry. Examples of application in- 
clude pipeline and airduct design as well as navigational raster 
charts. Although the paper is focused on the navigational use 
of raster grids, its results may be easily adapted elsewhere. 

The technology of electronic navigational charts (ENC) and 
digital electronic maps (DEM) has been used increasingly in 
maritime navigation, GPS and GIS systems [4]. As the raster 
charts are cheaper to be produced than vector ones they cover 
the majority of the electronic charts currently in use. In the 
ENC system one of the most common application is to deter- 
mine an optimal route between a start cell and a destination 
cell, that does not cross any obstacles (landmass, barriers, sho- 
als). The algorithm performing this should fulfill certain con- 
ditions. Of these the most obvious is that it does indeed find an 
optimal route if only such exists. Other conditions are those of 
needed time and memory space. In reality only algorithms of 
linear time and space complexities are useful as only such may 
be accepted by a real-time system processing large numbers of 
cells. The big O notation is further used for complexities, with 
O(n) indicating linear complexity [2]. 

The first solution to meet the above mentioned conditions 
was the maze routing algorithm presented by Lee [3], often 
described as a wave propagation process. To date, Lee’s algo- 
rithm and its variations are among the most widely used ro- 
uting methods finding applications in maze games, VLSI de- 
sign and road map routing problems. However, the original 
algorithm proposed by Lee has one serious drawback: it works 
only for the 2-geometry grid plane (also known as the Manhat- 
tan geometry). Only recently it has been upgraded to higher 
geometries while sustaining the linear time and space comple- 
xities. This new solution has been proposed by Chang, Jan and 
Parberry [1]. 


Despite the major progress, the potential use of the impro- 
ved Lee’s algorithm has still some limitations. Both the origi- 
nal algorithm and its upgraded version tend to find the shortest 
path which is not always identical to the optimal one. In pre- 
sence of many obstacles the algorithms determine a route con- 
taining so many turning points and course alterations that no 
navigator would be able to follow it even being so advised by 
a real-time routing system. Furthermore over larger geogra- 
phic areas both distance and number of turns will contribute to 
the total time spent traversing a tour. Thus minimizing the num- 
ber of turns is a desirable objective. Although the Chang, Jan 
and Parberry algorithm may cover the aspect of varied terrain 
(among other improvements) its data structure, based on that 
of the original Lee’s method, makes it unable to include the 
cost of a turn in path length. 

There is a number of methods (invented mostly for VLSI 
and automation purposes) that cover the problems of minimum 
bend path and shortest minimum bend path. Unfortunately, 
determining the shortest minimum bend path is of no value for 
long distance sea routing. Instead, the objective is finding the 
shortest paths with bend penalizing. The bend penalizing issue 
has been also considered in a number of works but the me- 
thods presented there are either not sufficiently fast or not ap- 
plicable for higher geometries. 

The article proposes a solution to this problem. A new data 
structure has been so designed as to reflect the cost of all cour- 
se alterations in each cell’s arrival time. An algorithm utilizing 
such structure has been implemented. The algorithm uses the 
following input parameters: user specified values of the costs 
of time of course alterations (turn penalties). It returns the de- 
termined path. 

In Sec. 2 the original Lee’s method and its version upgra- 
ded by Chang, Jan and Parberry are described in detail. Sec. 3 
presents a new algorithm that overcomes their limitations. In 
Sec. 4 some possible applications of this new solution are dis- 
cussed, and finally the conclusions are presented. 
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2. THE MAZE ROUTING ALGORITHM : 
PAST SOLUTIONS 


To facilitate checking out the Chang, Jan and Parberry ver- 
sion of the routing method its original notation is used here. 
In Sec. 3 the same notation is also used for the description of 
the new solution. 


Lee’s algorithm for the 2-geometry grid plane 


The popularity of this method lies in its simplicity and the 
guarantee to find the shortest path if only such exists. Since it 
is widely known no formal description is given here. 


Data structure 


The cell map static data is stored 
in an array containing two fields for every cell : 


SL — Boolean field whose value indicates whether a cell is as- 
signed to sea or landmass (or any other static obstacle). 
True for sea, and false otherwise. 

AT — Integer number field whose value is this cell arrival time, 
i.e. time needed to travel from the source cell to this cell. 


The dynamic data are stored in two lists : L; and L2 
containing cell pointers or indexes in the cell array. 


Algorithm overview 


A cell array is initialized with appropriate SL and AT valu- 
es. Two lists L; and L, are defined to keep track of the cells on 
the front of the wave (frontier cells) and their equal-distance- 
-step neighbouring cells, respectively. Putting the source cell in 
the list L; initializes the search. After all neighbouring cells of 
L; are included into the list L2, the list L2 is processed, so that 
an expanded wave front is found. Then every cell of L; is dele- 
ted if all of its neighbouring cells are processed (updated) and 
L; is updated by this new front of the wave. The search is ter- 
minated when the destination cell is found or — in the case of 
the barrier of obstacles — there are no more new cells to be 
processed, i.e. the last front of the wave has not generated any 
neighbouring obstacle-free cells. 


The Lee’s algorithm version 
upgraded by Chang, Jan and Parberry 


This method expands 2-geometry routing to any higher 
geometry. Details for an example 4-geometry routing are pre- 
sented beneath. Varied terrain aspect is not considered here as 
it increases the number of necessary lists, though it is generally 
supported by this version. 


2-geometry and 4-geometry movement rules 


Possible moves from a cell on the 2-geometry 
and 4-geometry grid planes are shown on the diagrams below. 


Fig.1. A 2-geometry neighbourhood 


Fig.2. A 4-geometry neighbourhood 


Data structure 
The static data is stored in an array containing 
three fields for every cell : 


SL — Integer number field whose value indicates whether a cell 
is assigned to sea or landmass (or any other static obstac- 
le). Its value is 1 for sea, and infinity for a landmass. 
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AT — Floating point number field whose value is this cell ar- 
rival time — time needed to travel from the source cell 
to this cell. 

VIS — Boolean field whose value indicates whether a cell has 
been already visited (inserted to temp-list) or not. True 
for a visited cell, false otherwise. 


The dynamic data is stored in three circularly used lists: L7, 
L2, L3 and one temporary list (temp-list). The index of a list 
currently in use is stored in a bucket_index variable. All lists 
contain cell pointers or indexes in the cell array. 


Algorithm overview 


There are some major differences between the 2 - and 
4-geometry algorithms. Firstly, additional lists are introduced 
to control the propagation speed in different directions (the ad- 
ditional two diagonal directions have single-step distance equ- 
al to /2 instead of 1). 

Secondly, the VZS cell field is introduced to make sure that 
each cell is inserted into temporary list exactly once. These 
two aspects make the algorithm keeping the original linear 
complexities. 


Algorithm formal description 


> The INSERT procedure inserts the index of c; into the temp- 
-list. 

> The procedure CLEAR empties any given list. 

> The RETRACE procedure retraces from the destination cell 
to the source cell according to their smallest AT value and 
builds a path ZL,,,, which is the output of the 4-geometry 
maze router. 


Algorithm (in Pascal) of 4-GEOMETRY-ROUTER 
(Cell-map, S, D, LL jain) 


Input: Cell-map, S, D 
Output: LL an 
begin 
bucket_index :=0; 
Loucket_index :=S ; 
VIS,: = TRUE ; 
temp-list : = @ ; 
path-exists : = FALSE ; 
while (Lpucket_index # Ø) or (Lbucket_index+1 # Ø) do 
if (D cell in Lpucket_index) then 


begin 
path-exists : = TRUE ; 
break while ; 
end 
for each cell c; in Lpucket index dO 
begin 
for each cell c; neighbouring c; do 
begin 
if (SL; = 1) then 
begin 
if (VIS, = FALSE) then 
begin 
VIS; : = TRUE ; 
INSERT(c;, temp-list) ; 
end 
Case |: 2-geometry neighbours 
AT pew = AT; +1; 
Case 2: diagonal neighbours 
AT ew : =AT; + J2 > 
if (AT new < AT;) then AT; (= AT new > 
end 
end 
end 


CLEAR (Lpucker index) 
if (temp-list + Ø) then 
begin 
for each cell c; in temp-list do 
INSERT(c; > Lfloor (Atj) mod 3) 
CLEAR (temp-list) ; 
end 
else bucket_index : = (bucket_index + 1) mod 3 ; 
end while ; 
if (path-exists = TRUE) then 
RETRACE(Cell-map(ATp) , LL pain) 
else path does not exist ; 
end . 


3. PRESENTATION OF THE NEW SOLUTION 


To simplify the presentation only the 4-geometry version is 
described here. However in general the algorithm is applicable 
to all geometries and it may also support varied terrain routing 
just the way the Chang, Jan and Parberry version does. 


Data structure 


The static data is stored in an array 
containing three fields for every cell : 


SL — Integer number field whose value indicates whether 
a cell is assigned to sea or landmass (or any other static 
obstacle). Its value is 1 for sea, and infinity for a land- 
mass. 

GAT- A sub-array of the floating point numbers. The size of 
the array is equal to the maximum number of the neigh- 
bouring cells and is 8 for 4-geometry. Each field of this 
sub-array contains the gate arrival times for different 
incoming gates of the current cell. Gate arrival time is 
the time taken to travel from the source cell to the cur- 
rent cell via certain neighbour of the current cell. 

VIS — Boolean field whose value indicates whether a cell has 
already been visited (inserted to temp-list) or not. True 
for a visited cell, false otherwise. 


The dynamic data is stored in circularly used lists: Z;...Ly , 
and one extra temporary list temp-list. The number n of the 
used lists strictly depends on the specified maximum course 
alteration cost and thus it is indirectly configured by means of 
algorithm input parameters. n = ceiling(maximum{single-step 
distances} + maximum{specified course alteration costs} + 1), 
where : maximum(single-step distances) is V2 for 4-geometry. 


Algorithm overview 


The key difference between the new solution and that pro- 
posed by Chang, Jan and Parberry is a replacement of each cell’s 
AT field with GAT array. The idea of incoming gates arrival 
times makes it possible to take into account course alteration 
costs without sacrifying the linear complexities that characteri- 
zed the previous versions. Every time cell data is updated, the 
arrival time of the appropriate gate is modified depending on 
the direction of the neighbouring cell that has initiated the update 
operation. The new candidate value of the gate arrival time 
field is determined according to the following formula : 


GAThew, j, gate_number — minimum {GATj,; + 
+ distance į ,j + delaygate number,1,--- GATi,g + 
+ distance į, jt delaY gate number,8}> 


where : 


i and j are indexes of the neighboring cells, 
gate_number is the current gate of the c; cell , 
numbers from | to 8 denote all gates of the c; cell. 


Delay values (penalties) are equal to zero for two gates of 
the same direction and have appropriate parameter values: d7, 
do or d3 for two gates whose direction difference is 45°, 90° or 
135°, respectively. The present GAT value is replaced with the 
candidate value if the new value is lesser than the current one. 


The below presented diagrams illustrate the way 
the GAT array values of the wave front cells are updated. 


Fig.4. A 45° turn 


Algorithm formal description 


* The procedures INSERT and CLEAR perform the same 
functions as those in the Chang, Jan and Parberry method. 

* The function GET_GATE_NUMBER (c;,c;) returns the 
number of the incoming gate of the cell c, , through which 
the travel from the cell c; is performed. 

* The procedure RETRACE retraces from the destination cell 
to the source cell and forms the output LL pan list. However, 
the rule of choosing the back way cells is different from 
that of the Chang, Jan and Parberry procedure. Here is cho- 
sen the cell whose sum of d; (i € {1,2,3}) modifier and the 
arrival time value of previous (closer to destination) cell’s 
incoming gate is minimum. 


Algorithm (in Pascal) of 44GEOMETRY-ROUTER-WITH- 
-TURN-PENALTIES (Cell-map, S, D, dı, dz, dz, LL yan) 


Input: Cell-map, S, D, dı 7 dz > d3 
Output : LLpath 
begin 
bucket_index:=0 ; 
L bucket index :=5 ; 
VIS; := TRUE; 
temp-list : =Ø ; 
path-exists : = FALSE ; 
all-lists-empty : = FALSE ; 
number-of-lists : = ceiling(./2+ maximum{d,, dz , d3}+1); 
while (all-lists-empty = FALSE) do 
if (D cell in Loucket_index) then 


begin 

path-exists : = TRUE ; 

break while ; 
end 

for each cell c; in Lpucket index dO 

begin 

for each cell c; neighbouring c; do 

begin 
if (SL; = 1) then 
begin 
if (VIS; = FALSE) then 
begin 
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VIS; : = TRUE ; 
INSERT(c,, temp-list) ; 
end 


Case 1: 2-geometry neighbours 
distance := 1 ; 
Case 2 : diagonal neighbours 
distance : = y2 ; 
gate_number : = GET_GATE_NUMBER(c;,c;) H 
GAT rew := GAT;, gate_number T distance > 
for each gate g, of the cell c; do 
begin 
Case 1: gate_number is 
the same direction as g% 
delay,:=0; 
Case 2: gate_number and 
g, difference is 45 degrees 
delay, : = dj); 
Case 3 : gate_number and 
g difference is 90 degrees 
delay, : = dz; 
Case 4 : gate_number and 
gx difference is 135 degrees 
delay, : = d3 ; 
GAT rey, k : = GAT;, k + 
+ distance + delay, ; 
if (GAT new, k< GAT new) then 
GAT new z= GAT pew, k> 


end 
if (GAT new < GAT;, gate mienber) then 
GAT;, gate_number SS GAT new > 
end 
end 
end 


CLEAR QLpneker index) 
if (temp-list + Ø) then 
begin 
for each cell c; in temp-list do 
begin 
k : = floor (minimum{ GAT; ; ,... GAT; , 3}) mod 
number-of-lists ; 
INSERT (c;, Lẹ) ; 
end 
CLEAR (temp-list) ; 
end 
else bucket_index : = (bucket_index + 1) mod 
number-of-lists ; 
all-lists-empty : = TRUE ; 
for each list L; do 
begin 
if (L; # Ø) then 
begin 
all-lists-empty : = FALSE ; 
break for 
end 
end 
end while ; 
if (path-exists : = TRUE) then 
RETRACE(Cell-map(GATp) , LLpatn) 
else path does not exist ; 
end . 


Computational complexity 


For each cell c; whose distance from the source cell is equ- 
al to or lesser than that of the destination cell, the following 
actions are performed : 


s% 


œ each of its neighbours c; is checked and possibly updated 


% for each of its neighbours c;, each of the gates of the cell c; 
is checked. 
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This gives the total of (2-A)-(2A)-n steps for the worst case, 
where A is a constant denoting geometry level (eight neighbo- 
urs and eight gates for 4-geometry) and n is the number of cells 
whose distance from the source cell is equal to or lesser than 
that of the destination cell. Thus the computational complexity 
of the proposed solution is O(n). 


Example results of application 
of the proposed routing method 
and the Chang, Jan and Parberry algorithm 


Below example results for the new routing method and the 
Chang, Jan and Parberry algorithm are presented. It has been 
assumed that costs of direction changes are: 1, 2 and 3 for 45°, 
90° and 135° degree turns, respectively. 


Fig.6. Route determined by using the proposed algorithm 
with direction change costs (penalties) dj = 1, dz = 2, d3 = 3 


In the above presented figures two different solutions of 
the same task are compared. The proposed algorithm (Fig.6), 
found a route with two turning points, different from that yiel- 
ded by the Chang, Jan and Parberry method (Fig.5) which de- 
termined a route with six turning points. Thus the total path 
lengths and penalties for both methods are : 


For the route by the Chang, Jan and Parberry method : 
basic path length = 8 + 5-/2 = 15.07 
total penalties = 5-d, + Il-dp =7 
total path cost = basic path length + total penalties =~ 22.07 


For the route by the proposed method : 


basic path length = 10 + 44/2 = 15.65 
total penalties = 2-d, = 2 
total path cost = basic path length + total penalties ~ 17.65 


From this comparison it can be stated that taking a seemin- 
gly longer (by 3.8%) route results in a much lower (by 22%) 
overall path cost. 


4. EXAMPLES OF APPLICATION 


The algorithm described in the paper is of general use. The- 
refore, whenever routing on the grid is involved and costs of 
direction changes matter, the benefits of using the proposed 
solution may be considerable. 

In navigation, real-time routing and collision avoidance 
systems are the most obvious fields where the algorithm may 
be adapted. The method is applicable for ENC systems and 
ECDIS display on raster plane in particular. For a detailed de- 
scription of a routing with collision avoidance system that may 
take advantage of the presented algorithm, the reader is refer- 
red to the paper [1]. 

Other potential maritime applications include pipeline and 
airduct routing, where minimizing the number of bends is cru- 
cial for cutting down the costs. However the algorithm would 
have to be first upgraded to a 3-D version. Fields of use outsi- 
de maritime industry may include VLSI design. 


5. SUMMARY 


In the paper a general searching method on the raster plane 
was presented. Owing to its new data structure the algorithm is 
capable of including costs of direction changes in the total cost 
of an optimal path, while keeping the linear computational time 
and space complexities. Therefore when costs of direction chan- 
ges are significant the solution is superior to those previously 
known. Possible applications include finding optimal ship rou- 
tes since their course alterations might be expensive or even 
risky in presence of obstacles. 


NOMENCLATURE 
List of variables used in the algorithms : 


all-lists-empty — indicates whether all buckets (lists) are empty 
or not 

— the arrival time (a time when the cell of the given 
index can be reached) 

— index of the bucket currently in use 


AT index 


bucket_index 


Cell-map — a map of cells, containing all relevant information 

d;,d,, ds — the turn penalties (delay times for 45, 90 and 
135-degree turns, respectively) 

delay — a delay resulting from the direction change 

distance — the distance between two neighbouring cells 


D — a destination cell 

gate number -— the number of the incoming gate, through which 
the cell is reached 

GAT ell index, gate number — the gate arrival time (a time when the 
incoming gate of the given cell can be reached) 

LL path — the returned path 

number-of-lists— the number of the buckets (implemented as lists) 

path-exists — indicates whether the searched path has been 
found or not 

S — a start cell 


SLindex — indicates whether the cell of the given index 
is sea (passable) or land (impassable) 

temp-list — a temporary list 

VIS — indicates whether a cell of the given index has 
already been visited 

Ø — empty list 


Abbreviations 


DEM -— Digital Electronic Maps 
ECDIS — Electronic Chart Display 
and Information System 


ENC — Electronic Navigational Charts 
GIS — Geographical Information System 
GPS — Global Positioning System 

VLSI — Very Large Scale Integration 


2- and 4-geometry — two, four trajectories, 
respectively 


3-D — three dimensions 


BIBLIOGRAPHY 


1. Chang K. Y., Jan G.E, Parberry I. : A Method for Searching 
Optimal Routes with Collision Avoidance on Raster Charts. 
Royal Institute of Navigation. The Journal of Navigation 
No.56/2003 

2. Kubale M. : Introduction to Computational Complexity. 
Wydawnictwo Politechniki Gdanskiej (Publishing House of 
Gdansk University of Technology). 1994 

3. Lee C.Y. : An Algorithm for Path Connection and Its 
Applications. IEEE (Institute of Electrical and Electronics 
Engineers) Transactions on Electronic Computers, EC-10, 1961 

4. Weintrit A. : The Electronic Chart Systems and Their 
Classification. Annual of Navigation, No 3/2001. Polish 
Academy of Sciences. Polish Navigation Forum, Gdynia 2001 


CONTACT WITH THE AUTHOR 


Rafat Sztapezynski, M.Sc.,Eng. 
Faculty of Ocean Engineering 
and Ship Technology, 
Gdansk University of Technology 
Narutowicza 11/12 
80-952 Gdansk, POLAND 
e-mail : rafal@pg.gda.pl 


aa C onference — 
SCIENTIFIC SEMINAR OF 
REGIONAL GROUP 
of the Section 
on Exploitation Foundations 


On 25 November 2004 the Institute of Fluid-Flow 
Machinery, Polish Academy of Sciences (PAS), Gdańsk, 
hosted the last-in-the-year scientific seminar of the Re- 
gional Group of the Section on Exploitation Foundations, 
Machine Building Committee, PAS. 

During the main part of the seminar 4 papers were pre- 
sented whose authors came from scientific staff of the In- 
stitute, namely : 


+ On energy and entropy — by J. Mikielewicz 

+ Plasma techniques applied to modern cars 
by J. Mizeraczyk 

+ Experimental modal analysis — by M. Łuczak 

+ Structural funds — research projects and 
their financial assembling — by J. Kiciński 


After interesting discussion on the papers Prof. J.Gir- 
tler, the Group's Chairman, presented a proposed plan of 
seminars for the year 2005. 

The seminar was ended by presentation 
io the up-to-date laboratory facilities of the Institute. e 
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OPERATION & ECONOMY 


FOREIGN 


a 


PSAM 7 — ESREL’04 


On 14-18 June 2004 in Berlin had place 
the International Conference on : 


Probabilistic Safety Assessment 
and Management 


It was organized by two bodies : International Associa- 
tion for Probabilistic Safety Assessment and Management, 
and European Safety and Reliability Association. 

Due to such arrangement its participants took part in 
a very broad meeeting of scientists, researchers and experts 
of many branches, which was demonstrated in 594 papers 
prepared by authors from 40 countries, from all the world 
around. 

Among the gremium were also representatives of Polish 
scientific research centres who presented the following to- 
pics : 
= Environmental safety of sea-going ship power plant — by 

Brandowski A., and Liberacki R. (Gdansk University of 

Technology) 
= Risk in the aspect of safety and reliability of autonomous 

system — by Drobiszewski J. and Smalko Z. (Warsaw 

C University of Technology) 


=> New methods of solving multidimensional non-linear 
optimization problems — application thereof to construct 
air-line flight schedules — by Jazwinski J., Klimaszew- 
ski S., and Zurek Z. (Air Force Institute of Technology, 
Warsaw) 

= Availability improvement of port transportation structu- 
res and their operation processes — by Kotowrocki K. 
(Gdynia Maritime University) 

= Incorporation of human and organizational factors into 
qualitative and quantitative risk analyses- by Kosmow- 
ski K.T. (Gdansk University of Technology) 

= Modelling and uncertainty in system analysis for safety 
assessment — by Kosmowski K.T. (Gdansk University 
of Technology) 

= On reliability improvement of the port grain transport 
systems — by Kwiatuszewska — Sarnecka B. (Gdynia 
Maritime University) 

© Reliability model of combined transportation system 
by Nowakowski T. (Wrocław University of Technology). 


Moreover, Mr J. Gorski of Gdansk University of Tech- 
nology took part in the activity of Application Area Co-ordi- 
nators, a 16-person international organizational team, within 
the area : Information Technology and ee en 


Miscellanea 


Without doubt welding is one 
of the most important fields of engineering. 


Welding, soldering and bonding — these are engineering 
processes which play an important role in production, repair 
and operation of objects. Therefore welding is teached in 
many engineering schools and universities worldwide. 


However it has turned in practice that their graduates — 
— welding engineers with their qualification do not quaran- 
tee to cope with required responsible supervision over reali- 
zed welding processes, which detrimentally influence quali- 
ty and durability of produced engineering objects. In order 
to improve the situation European Federation for Welding 
(EFW) and International Institute of Welding (ITW) came 
to the conclusion that welding processes — from the point of 
view of quality assurance — require a.o. special surveillance 
by specially trained engineering personnel. 


Education of such personnel is carried out during a spe- 
cial training supplementary to academic studies. It is reali- 
zed in accordance with the requirements defined in guide- 
lines for unified educational system (Minimum Requirements 


World qualification level ES 
è 
for welding engineers 


for Education, Examination and Qualification) issued by the 
above mentioned organizations. The requirements contain 
overall description of aims and scope of the education, as 
well as conditions and criteria for carrying out final exami- 
nation to obtain the diploma of European/International Wel- 
ding Engineer (EWE/IWE). 


In Poland it was Gdansk University of Technology 
which first undertook such a task, where at its 
Mechanical Faculty the IIW Approved Training Body 
was organized. 


Its activity was initiated during spring semester of 2001/ 
2002 academic year under the authority of Welding Institu- 
te, Gliwice, being the IW Authorized National Body. The 
first 300-hour course of the education on the highest inter- 
national level was completed in July 2004. All 24 partici- 
pants of the course passed the very difficult final exam with 
the mean grade of "almost very good". This way they ente- 
red the group of welding engineers of the highest internatio- 
nal level which makes it possible to take up the highest posts 
in this profession. 
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